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ABSTRACT 

We have measured the HI power spectrum of the nearly face-on spiral galaxy 
NGC 1058 from radio-interferometric observations using a visibility based es- 
timator. The power spectrum is well fitted by two different power laws P(U) = 
AU a , one with a = — 2.5±0.6 at small length-scales (600 pc to 1.5 kpc) and 
another with a = — 1.0 ± 0.2 at large length-scales (1.5 kpc to 10.0 kpc). We 
interpret this change in the slope of the power spectrum as a transition from 
3D turbulence at small length-scales to 2D turbulence in the plane of the 
galaxy's disk at large length-scales. We use the observed break in the power 
spectrum to estimate the galaxy's scale-height, which we find to be 490 ± 90 
pc. 
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1 INTRODUCTION 

Evidence has been mounting in recent years that turbulence plays an important role in deter- 
mining the physical conditions of the neu tral interstellar medium (ISM) as well as generating 
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the hierarchy of structures seen in it (see 
for recent reviews). Observational evidence includes the fact that the fluctuation power 
spectrum of a variety of tracers (HI 21-cm emission intensity, HI 21-cm optical depth, dust 
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and molecular emission) is a scale-free power-law (ICrovisier fe Dickey 
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This scale-free behavior is characteristic of a turbulent medium. Similarl y, the HI distri- 



butio n in several dwarf galaxies in the M81 group appears to be fractal (jWestpfahl et al. 



19991 ). Further, the Fourier transform power spectra of the V and Ha images of a sample of 



irregular galaxies are also found to be a power law, indicating that there is no characteristic 
mass or luminosity scale for OB associations and star complexes (Willett et al. 2005). 



Recently iBegum et al.l (120061 ) have presented a visibility based formalism for determining 
the power spectrum of HI intensity fluctuations in galaxies with extremely weak emission. 
This was applied to the dwarf galaxy DDO 210 and the power spectrum was found to be 
a power la w with slope —2.8 ± 0.4 over the length scales 100 — 500 pc. In a subsequent 



paper, we (IDutta et al.ll2008l ) have applied the same formalism to an external disk galaxy 



NGC 628 (M 74) and found the slope to be -1.7 ± 0.2 over the length scales 0.8 - 8.0 kpc. 
This is significantly shallower than the power spectra in earlier measurements including our 



Gal axy, the SMC and DDO 210, a 
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of w 
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lich exhibit a slop e roughly in the range —2.5 to 



Begum et al.ll2006l ). The shallower slope observed 



in NGC 628 was interpreted as arising from two dimensional (2D) turbulence in the plane 
of the galaxy's disk on length-scales larger than the galaxy's scale-height. The steeper slope 
seen in the other galaxies was ascribed to three dimensional (3D) turbulence at small length- 
scales. We expect a transition from 3D to 2D turbulence at a length-scale corresponding to 
the galaxy's scale height. An observational de tection of this trans i tion provides a method to 



estimate the scale-height of face-on galaxies (IPadoan et al. 



2001 



Elmegreen et al.ll200ll ). a 



quantity which is otherwise difficult to measure. 



NGC 1058 is an almost face-on, type Sc, spiral galaxy with an inclination angle in 
the range 4° to 11° ( IPetric fc Rupenl 120071 ) . In this Letter we use a break in the observed 
HI power spectrum to determine the galaxy's scale-height. The distance t o this galaxy 



i s uncertain with preyious e stimates ranging from 10.0 Mpc (IBorosonl 119811 ) to 14.5 Mpc 
( ISandage fc Tammannl Il974l ). Throughout this paper we adopt a distance of 10 Mpc for 
NGC 1058. At this distance 1" corresponds to 48.5 pc. 
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Figure 1. The 13.6" X 13.6" resolution integrated HI column density map of NGC 1058 (contours). The contours levels are 
3., 10., 30., 50., 70., 100., 120., 130., and 140. xlO 20 cm" 2 . 

2 DATA AND ANALYSIS 

We have used archival HI data of NGC 1058 from the Very Large Array (VLA). The obser- 
vations had been carried out on 14 th June 1993 in the C configurations of the VLA. The data 
was downloaded from the VLA archive and reduced in the usual way using standard tasks in 
classic AIPsE. The HI emission from NGC 1058 spans 24 channels (i.e, channel 53 to channel 
76) of the 127 channel spectral cube. We have used only the central 16 channels (i.e, channel 
57 to channel 72) with relatively strong HI emission. The frequency width of each channel 
corresponds to 2.58 km s _1 . The continuum from the galaxy was subtracted from the data 
in the uv plane using the AIPS task UVSUB. The resulting continuum subtracted data was 
used for the subsequent analysis. Figure [1] shows a total HI column density (Moment 0) map 
of NGC 1058 made from this data. The HI disk of the galaxy is nearly face-on. The angular 
extent of the HI distribution in Figured] is measured to be 1 1/ x 13' at a column d ensity of 
10 19 atoms cm -2 , which is ~ 4 times it's Holm berg diameter (IPetric fc Rupenl 120071 ) . 



Begum et al. 



f|2006h and 



Dutta et al 



(120091 ) contains a detailed discussion of the visibility 
based HI power spectrum estimator Pm(U) = ( V U (U)V*(U + AU) ), hence we present only 
a brief discussion here. Here U refers to a baseline, the antenna separation projected in 
the plane perpendicular to the direction of observation, measured in units of the observing 
wavelength A. Note that, U expressed this way is a dimensionless quantity. It is common 
practice to express U in units of kilo wavelength (kA). Every visibility V V (U) is correlated 
with all other visibilities V*(U + AU) within a disk \ AU\ < (ttOo)" 1 , where 8q is the angular 



NRAO Astrophysical Image Processing System, a commonly used software for radio data processing. 




U (kA) 

Figure 2. Real (blue, solid) and imaginary (green, dot-dash) parts of the observed value of the HI power spectrum estimator 
Pffl(U) f° r N = 16. lcr error-bars are shown only for the real part. The real part of the Pm(U) from line-free channels (red, 
bold dot-dash) is also shown. 



extent of the galaxy. The correlations are averaged over different U directions assuming 
that the signal is statistically isotropic in the plane of the galaxy's image . To increase the 
signal to noise ratio we further average the correlations in bins of U and over all frequency 
channels with HI emission. The expectation value of the estimator Pm{U) is real, and it is 
the convolution of the HI power spectrum Pbi(U) with a window function |W„({7)| 2 . The 
window function quantifies the effect of the large-scale HI distribution and the finite angular 
extent of the galaxy, and it can be assumed to be sharply peaked around U — with a width 
of order (ttOq)^ 1 . At baselines U ^> (ttOq)^ 1 the function iH^f/)! 2 is well approximated by 
a Dirac delta function and the expectation value of Pm(U) gives an estimate of the HI 
power spectrum Pm(U). The effect of |W^(i7)| 2 in estimating the power spectrum and the 
considerations for determining (tiOq)^ 1 is discussed in detail in iDutta et al.l (120091 ) . 

The 1 — a error-bars for the estimated power spectrum is a sum, in quadrature, of 
contributions from two sources of uncertainty. At small U the uncertainty is dominated by 



the sample variance which comes from the fact that we have a finite and limited number 
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of independent estimates of the true power spectrum. At large U, it is dominated by the 
system noise in each visibility. 

To determine if the slope of the HI power spectrum changes with the width of the 
frequency channel, we have combined N successive channels to obtain a data set with 16 /N 
channels of width A2.58 x km s _1 each. We have determined the HI power spectrum for a 
range of iV values. 



3 RESULTS AND DISCUSSION 



Figure [2] shows the real and imaginary parts of PsiiU), which is th e observed value o f 
the estimator Pm{U). As expected from the theoretical considerations (IBegum et al.ll2006l ). 
the imaginary part is well suppressed compared to the real part. To test for a possible 
contribution from residual continuum, we also show the real part of Pri(U) estimated using 
line free channels. This is found to be much smaller than the signal. For the channels with 
HI emission, the observed Pm{U) may be directly interpreted as the HI power spectrum at 
U values that are considerably larger than (tiOq)^ 1 = 0.1 kA. 

The power spectrum is well fitted by two different power laws P(U) = AU a , one with 
a = — 1.0 ± 0.2 for U = 1.0 kA to 6.5 kA (large length-scales) and another with a = 
— 2.5 ± 0.6 for U = 6.5 kA to 16.0 kA (small length scales). The results are tabulated in 
Table [TJ The best-fit power-law were dete rmined by a y 2 mi nimization with respect to the 
amplitude, A and the power-law index, a (IBegum et al.ll2006l ). The I — a error-bar on a, the 
parameter of interest here, is determined by projecting the ell ipse corresponding to A% 2 = 1 



Press et al. 



Jl992J)). 



onto the a axis in the a — A parameter space (page 694 of 

Both HI density fluctuations as well as spatial fluctuations in the velocity of the HI 
gas contribute to fluctuatio ns in the HI specific intensity. Considering a turbulent ISM, 
Lazarian fc Pogosyan (J2OOO) have shown that it is possible to disentangle these two contri- 
butions by studying the behavior of the HI power spectrum as the thickness of the frequency 
channel is varied. The slope of the HI power spectrum is expected to change with t he fre - 
quency channel th ickne ss if it is du e to th e gas velocities, However, IBegum et al.l (120061 ). 



Dutta et al. 



(120081 ) and iDutta et al.l (120091 ) do not fi nd this change fo r any of the galaxies 



in their sample of spiral and dwarf galaxies. Further, 



Elme green et al 



(120011 ) also reported 



a similar behavior for LMC We find that the HI power spectrum of NGC 1058 does not 
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Baseline Range 


N 


AV 


a 


(kA) 




(km s" 1 ) 






1 


2.68 


-1.0 ±0.2 


1.0 - 6.5 


8 


21.44 


-0.9 ±0.2 




16 


42.88 


-0.8 ±0.2 




1 


2.68 


-2.5 ±0.6 


6.5 - 16.0 


8 


21.44 


-2.2 ±0.6 




16 


42.88 


-2.2 ±0.5 



Table 1. NGC 1058 power spectrum are well fitted by two power laws of Phi = AU a in two different baseline ranges. This 
table summaries the result for different channel width. 



exhibit a statistically significant change with increase of the channel thickness (Table [I]), 
indicating that the observed power spectrum is due to HI density fluctuations. 

A baseline U corresponds to an angular scale of 1/U in the plane of the sky. In the 
subsequent discussion, we use D — 10 Mpc (distance to the galaxy) to convert the baseline 
U to a length-scale D/U in the plane of the galaxy's image. We find that the slope of the 
power spectrum of NGC 1058 is a = — 2.5±0.6 at small length-scales (0.6 — 1.5 kpc) whereas 
it is a = —1.0 ±0.2 at large length-scales ( 1.5 — 10.0 kpc). The largest lengt h-scale (10 kpc) 



is definitely larger than the t ypical HI scale heights within the Mi 



198 



4; 



Wouterloot et al 



19901 ) and in external spiral galaxies (e.g. 



ky-Way (Lock man et al. 



Narayan fc Jog 



2002 ). It 



is then quite reasonable to conclude that the slope a — — 1.0 ± 0.2 is that of 2D turbulence 
in the plane of the g alaxy's disk. Simu lations of the power spectrum of HI emission from a 
face-on disk galaxy (jDutta et al.ll2009i ) predict a change in the slope a corresponding to a 
transition from 2D to 3D turbulence. This is expected to occur at a baseline U = D/irzh, 
where is the galaxy's scale-height. Interpreting the change in the slope of the observed 
power spectrum seen at U = 6.5 kA as the transition from 2D to 3D turbulence, we find the 
scale-heigth of NGC 1058 to be 490 pc. Further, we interpret the slope of a = — 2.5 ± 0.6 to 
be that of 3D turbulence. To our knowledge this is the first observational determination of 
the scale-height of a nearly face on spiral galaxy through its HI power spectrum. The width 
of the bin AU, corresponding to the baseline where the break occurs, gives an estimate of the 
uncertainty in the e stimated scale- height as Azh = DAU /ttU 2 , whereby we have 490 ±90 pc. 
Kregel et al.l (12004 ) present HI images of a large sample of edge on intermediate to late type 
spirals; from their data the ratio of the HI disk height to the radius of the HI disk (at a 
column density of lM pc -2 ) is ~ 0.06 ± 0.015. Since the disk of NGC 1058 has a diameter 
of 35 kpc (Figured]) at this column density, we expect it's scale-height to be 1.05 ±0.26 kpc, 
in excess of our result. Note that the HI scale-height is known to increase with the distance 
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from the center, and our result is an average value over the entire disk with possibly a larger 
contribution from the inner parts where the emission is strongest. 

A similar behavior of a change in the slope of the power spectrum at larger length scales 
has been also observed in LMC (Elmegreen et al. 2001). The HI power spectrum of the LMC 
flattens at large length-scales, which was interpreted as a transition from three-dimensional 



to two- dimensional turbulence . Using this method t 



100 pc (jElmegreen et ai- 



re scale-height of LMC w as found to be 



200 ll ). It is well known (eg. 



Elmegreen fc Scalo 



200J) that the slope 



of the velocity power spectrum changes by unit y from —11/3 to — 8/3 in going from 3D to 
2D for incompressible Kolmogorov turbulence (jKolmogorovl 1 19411 ) . This may be related to 
the change in slope from —2.5 to —1.0 measured in NGC 1058, though the turbulence here 
is compressible and the observations are those of density fluctuations and not the velocity. 
The relation between the power spectrum of density fluctuations and the velocity power 
spectrum is not very clear for compressible turbulence. 

There is now mounting evidence that spiral galaxies exhibit scale-invariant density fluc- 
tuations that extend to length-scales of ~ lOkpc (eg. NGC 628, NGC 1058) which is com- 
parable to the diameter of the HI disk. While a large variety of possible energy sources 
l ike proto-stellar winds, s upernovae, shocks, etc. have been proposed to drive turbulence 
(jElmegreen &: Scald 120041 ). it is still to be seen whether these are effective on length-scales 
as large as 10 kpc. 

It is interesting to note that a break in the power spectrum, like the one seen here 
for NGC 1058, h as also been observed in the shell-type supernovae remnant Cassiopeia A 



teoy et al. 



20081 ). The slope of the power spectrum changes from —2.22 ± 0.03 at 1.6 — 10 kA 



(large length-scales) to —3.23 ± 0.09 at 11 — 30 kA (short length-scales). This change in the 
slope is interpreted as a transition from 2D to 3D magneto hydrodynamic turbulence. The 
transition occurs at a length-scale of 50 pc which corresponds to the thickness of the shell. 
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